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Abstract 

We describe the spectrum of a non-self-adjoint elliptic system on a finite 
interval. Under certain conditions we find that the eigenvalues form a dis- 
crete set and converge asymptotically at infinity to one of several straight 
lines. The eigenfunctions need not generate a basis of the relevant Hilbert 
space, and the larger eigenvalues are extremely sensitive to small pertur- 
bations of the operator. We show that the leading term in the spectral 
asymptotics is closely related to a certain convex polygon, and that the 
spectrum does not determine the operator up to similarity. Two elliptic 
systems which only differ in their boundary conditions may have entirely 
different spectral asymptotics. While our study makes no claim to general- 
ity, the results obtained will have to be incorporated into any future general 
theory. 
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1 Introduction 



Consider a first order elliptic differential system acting in L 2 := L 2 ((a, (3), C n ) 
according to the formula 

(Lf)(x):=A(x)f(x) 

where A is an n x n complex matrix depending on x. (We show in Section 2 
that our methods also apply to certain second and higher order elliptic systems.) 
If A(x) is a piecewise continuous function of x then the solution to the equation 
Lf = zf can be written in the form 

f(x) = U(z,x)f(a) 

where U(z, x) is an invertible n x n matrix which depends continuously on x and 
analytically on z. The operator has domain contained in W 1,2 := W 1,2 ((a, /?), C n ). 
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Since W 1,2 C C[a,/3], we can impose general boundary conditions of the type 

Sf(a) + Tf((3) = 0. (1) 
where S, T are any n x n matrices. 

Problems of this type arise in several contexts. One relates to non-equilibrium 



thermodynamics, and another concerns the study of generalized determinants 



for elliptic operators on manifolds with boundary, |I3[ [T7|, [18]]. Of course problems 
involving higher order differential operators can also be reformulated in terms of 
first order systems; we do not, however, assume Hamiltonian structure, which 
would force n to be even and also imply strong constraints on the solutions of the 
equations. The present paper provides a class of nearly exactly soluble examples 
which illustrate some of the phenomena such studies have to face. 

The spectral behaviour of L has obvious relevance to the time-dependent system 

If each matrix A(x) has only real eigenvalues a r (x) then (|2]) is called a hyperbolic 
system and a r (x) are the characteristic speeds at x. This case has special features 
by comparison with the case in which the eigenvalues are in general position in 
the complex plane, and is only partially analyzed in this paper. See Section 7 on 
semigroup properties. 

Example 1 Consider the 'Dirichlet' boundary conditions f{a) G U, f((3) G V, 
where U, V are linear subspaces in C n with 

dim({7) + dim(V) = n. 

This falls within the above scheme if we choose S, T as follows: S should have kernel 
U and range W while T should have kernel V and range X where W D X = {0} 
and W + X = C n . In particular if U is the linear span of the single vector u ^ 
and V is the annihilator of the vector d / then the boundary value problem 
reduces to 

F(z) := (U{z,(3)u,v) = 0. 

Returning to the general context, the embedding of W 1 ' 2 into L 2 is compact, so 
the operator L has compact resolvent if its spectrum is not equal to C, and its 
spectrum must consist of discrete eigenvalues of finite multiplicity. If A{x) = al 
for all x G [a,0\ then the only possible solution of Lf = zf is f(x) = ce zx ^ a 
for some c G C n ; in the context of Example [l] we deduce that Spec(L) = C if 
U fl V 7^ {0}, but in all other cases the spectrum is empty. The situation changes 
entirely if we assume that A(x) is constant but its eigenvalues a r are all different 
and all non-zero. 

It follows immediately from the above discussion that z is an eigenvalue of L if and 
only if 

F(z) := det(S + TU(z,0)) = 
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where F is an entire function of z. There are well-developed numerical procedures 
for evaluating F(z) for any z and for finding the points at which F vanishes. 
However, the asymptotic behaviour of the eigenvalues of L and the study of such 
quantities as the regularized determinant are difficult subjects. There is one case in 
which they can be analyzed fairly directly, namely when A(x) is piecewise constant, 
an assumption which we make for the remainder of the paper. 

In the following sections we present a detailed analysis of the function F(z). We 
prove in the 'generic case' that there are several series of eigenvalues diverging to 
infinity along straight lines which are determined by a certain convex polygon. The 
asymptotic form of the spectral counting function depends upon the length of the 
boundary of this polygon. There are exceptional cases in which almost periodic 
structure arises, and which we do not analyze fully. The implications of our results 
for the spectral analysis are explained with a series of simple examples. 

A much more complete analysis of the second order case with n = 2 has been 
performed by Boulton, [|, ||, and his results have informed our approach to the 
case n > 2. Our results describe how the more general analysis goes in the generic 
case, and what further phenomena have to be considered for a complete analysis. 

2 The function F(z) 

We assume henceforth that 

a = «o < di < ... < a m = (3 

and that A(x) = A s if a s -i < x < a s . We also assume that each A s is invertible 
and diagonalizable. Then 

U(z,/3) = U m U m - 1 ...U 2 U 1 

where 

and each matrix D s is diagonal. We then have 

F{z) = det(S + TV m e zD -V- 1 .V m .ie zD ^V m ^ l ...Vie zD ^V 1 - 1 ) 

= X>e*> (3) 

r=l 

It might be thought that the case in which A(x) is continuous can be obtained from 
the above by an approximation procedure. In Example |9| we show that taking such 
a limit would not be a straightforward matter. We discuss the spectral asymptotics 
of such piecewise constant operators in Section 4, and make a conjecture about the 
case in which A(x) depends continuously on x. 
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Let us look at the case of Dirichlet boundary conditions in more detail. Assuming 
that A is independent of x, invertible and diagonal with eigenvalues a r , the solution 
of the eigenvalue equation is 

f r {x) = c r e z{x - a)/ar 

where c G C" must be chosen so that / satisfies the boundary conditions. If 
u n - p is a basis of the annihilator of U in the dual space and Vi,...,v p is 
a basis of the annihilator of V then the conditions are that (c, = for all 
1 < i < n — p and (c, Wj) — for all 1 < j < p, where Wj G C n is defined by 

Wj ^ <p — ^ j jf ^ * 

The existence of a non-zero eigenvector then reduces to 

det{vi, v n - p , wi, w p } = 0. 

This may be rewritten in the form F(z) = where F is defined by (f|). In this 
case each of the coefficients 7 r is a sum of p distinct terms of the form (/3 — a) /aj. 

In the simplest case dim({7) = 1 and there are n distinct 7 r . In general the number 
of 7 r is given by a combinatorial expression. However, some of the 7 r may coincide, 
and some of the 5 r may vanish. 

We next consider a similar problem for second order elliptic systems. Let H act in 
L 2 := L 2 ((a,j3), C n ) according to the formula 

(Hf) r (x) := 

for 1 < r < n, where a r ^ for all r. The operator has domain contained 
in W 2 ' 2 := W 2 ' 2 ((a,/3),C n ). Since W 2 ' 2 C C«[a,/3], the boundary conditions 
may involve /(a), f'(a), /'(/?)• We do not consider the most general choice 

of boundary conditions, which leads to a somewhat complicated function F(z), 
but follow Boulton in assuming a combination of Dirichlet or Neumann boundary 
conditions in the following sense, [|, ||. We suppose that /(a) G Z7i, e ^2, 

G Vi and f'((3) G V 2 where U\, U 2 , Vi, V 2 are linear subspaces of C n . We 
impose the minimal further condition 

dim(L r i) + dim(C/ 2 ) + dim(14) + dim(V 2 ) = 2n. 

The eigenvalue problem may be written in the form Hf = z 2 f, and its solution is 
of the form 

f r (x) = c r e {x - a)z/ar + d r e- {x - a)z/ar 

where (c, d) G C 2n , provided z 7^ 0. The case z = is dealt with separately. The 
boundary conditions lead to 2n linear constraints on the vector (c, d) which have 
a non-zero solution if F(z) = 0, where F(z) is a function of the form (^). 

As a particular case we mention the possibility H := L 2 where L is the first order 
operator already discussed. This corresponds to the choice of boundary conditions 
f(a) G U, f'(a) G A~ l U , f(J3) G V, f'{(3) G A~ X V, where A id = 5^. The 
eigenvalues of H are the squares of the eigenvalues of L, and it follows from our 
analysis below that they are asymptotic to certain parabolae at infinity. 
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3 Zeros of F(z) 



Our task is to study the asymptotic distribution of the zeros of entire functions of 
the form 

F(z)=jT6 r e^ = j:F r (z) (4) 

r=l r=l 

as \z\ — > oo. We assume that all 7 r are different, that all 5 r are non-zero and 
that R > 2. We will see that the asymptotic structure of the set of zeros depends 
heavily on the convex hull K of the set of j r . We will only consider the generic 
case, defined as follows. After relabelling, the convex hull has vertices 7i,...,7q 
written in anticlockwise order, while the remainder of the 7 r lie in the interior of 
K; we assume that no j r lies on a edge of K unless it is already a vertex of K. If 
1 < r < Q we define 

if 1 < r < Q - 1 
if r = Q, 

if 2 < r < Q 
if r = 1. 

We will prove that, in addition to a finite number of zeros near to the origin, F(z) 
has Q series of zeros, each of which converges asymptotically at infinity towards 
one of the sets on which 

\F r (z)\ = \F r _(z)\. 
where 1 < r < Q. Rewriting this in the form 

l^le*^ = |5 r _|e^- 

we see that it is a straight line 

z ■ (7r - 7r_) = K (5) 
perpendicular to the edge (7 r _,7 r ) of K. 

If one or more of the 7 r with Q < r < R lie within an edge of K then the spectral 
behaviour of L at infinity is much more complicated. Boulton has shown that the 
behaviour depends upon whether those 7 r which lie within an edge of K divide 
that edge into parts whose lengths have rational or irrational ratios. In the former 
case the eigenvalues still converge at infinity to one of several straight lines, but 
there are several lines perpendicular to the relevant edge instead of just one. In 
the irrational case the eigenvalues have almost periodic structure at infinity within 
one of several strips, with one strip perpendicular to each edge. Boulton's analysis 
is only presented in a typical second order case, but the ideas clearly extend as 
stated, [|, I . 

Returning to the generic case, we study the asymptotic location of the zeros of 
F(z) by dividing the exterior of K into 4Q subregions. Given e > there are 3 



r + 1 
1 

r — 1 

Q 
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semi-infinite strips S^, S+, S~ associated with each edge (7 r _,7r) of K, and also 
a wedge W r associated with each vertex r ) r . Let e r be the outward pointing unit 
normal to the edge (7 r _ , 7 r ) . Given 1 < r < Q we define 

SI := {z G C : z • e r > 7 r • e r and k r — e < z ■ (7,, — j r _) < k r + e} 

:= {z G C : z ■ e r > 7 r • e r and k r + e < z • (7,, — 7 r _) < k r + 1} 

57 := {z 6 C : z ■ e r > 7, • e r and fc r — 1 < 2; • (7 r — 7 r _) < A; r — e} 

W r := {z K : k r + 1 < z ■ (j r — 7 r _) and k r+ — 1 > z • ( r y r+ — 7 r )}- 

Theorem 2 Given < e < 1/9 i/iere exists T £ such that any eigenvalue z of L 
satisfying \z\ > T £ must lie in S^. for some r such that 1 < r < Q . 



If s 7^ r, r_ then (7,. — 7 S ) • e r > and for large enough \z\ we have 



A similar argument applies to S r . 
If s ^ r, r± then either 

lr-ls= P{lr - 7r_) + Q^r 

where p > and g > 0, or 

7r - 7s = P(7r - 7r+) + qe r+ 
where p > and g > 0. It follows that 



Proof If z G 5+ then 



F r _(*)/F r (z)| <e" e < l-e/2. 



F a (z)/F r (*)| <e/2i2. 



Combining these bounds we see that for such z 



F(z)\ > \F r (z)\e(l/2 



(R-2)/2R) > 0. 




m „. {{lr-ls)-z} = 

,Z€W r 



Now let z eW r . We have 



F r± (z)/F r (z)\ <e 



-1 



and if 1 2; I is large enough we also have 



F s (z)/F r (z)\ <e/2R. 



provided s 7^ r,r±. Combining these inequalities we obtain 



F(z)\ > \F r (z)\(l-2e 



-1 



e(R-3)/2R) > 0. 
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Note 3 If one of the 7 S for Q < s < R is very close to the boundary of K then 
T £ will be correspondingly large in the above proof. This is inevitable because the 
possible asymptotic directions of the spectrum change as 7 S moves through an edge 
of K, whereupon K changes so that 7 S becomes another vertex of K. 

Note 4 The semi- infinite lines (||) divide the complex plane into Q sectors S r . If 
\z\ — > oo within S r then F(z) ~ S r e z/yr , so 

\og(F(z)) ~ log(5 r ) + 277. 

This formula may be used to compute the regularized determinant in the sense of 
PI, ITTL pi. 



4 Spectral Asymptotics 

We have proved that the zeros of F(z) converge asymptotically towards one of Q 
straight lines as \z\ — > oo. To obtain more detailed information we introduce the 
functions 

G r {z) := F r {z) + F r _{z). 

Theorem 5 The zeros of G r (z) lie on the line ([Jj with constant distance 2n/\ / y r — 
7 r _ | between any two consecutive zeros. 

Proof If we put 7 r — ■y r _ = pe 10 where p > and 6 e R and c r = log(S r _/5 r ), then 
the zeros of G r (z) are given by 

z = (cv + (2n + l)Tci)e ie /p 

where nfZ. The statements of the lemma all follow from this. 

Theorem 6 The zeros of F(z) which lie in converge as \z\ — > oo to the zeros 
of G r (z), in the sense that the modulus of the differences of corresponding zeros 
converges to zero. 

Proof We have already proved that the zeros of F(z) converge to one of the lines 
(|5|). The statement is a straightforward application of Rouche's theorem, since 
we have already noted that the remaining terms in the series are asymptotically 
negligible as \z\ — » oo within S^. 

Corollary 7 If N(E) is the number of zeros of F(z) such that \z\ < E then 

N(E) = b(K)E/2n + 0(l) 
as E —>■ oo, where b(K) is the length of the boundary of K. 
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Proof The number of zeros of F(z) satisfying \z\ < E associated with each line is 
\lr — It- \E/2tt + 0(1) by the above theorem, since the zeros of F(z) only converge 
to the zeros of G r (z) in one direction. The corollary follows by summing over r. 

Example 8 If Lj are two elliptic systems acting in L 2 ((a, 0), C ni ) respectively for 
i = l,2 then we may consider their direct sum L = Li ®L 2 . By applying Corollary 
[7| to each of the components we obtain 

N(E) = (6(#i) + b(K 2 ))E/2ix + 0(1) 

in an obvious notation. On the other hand we have F(z) = Fi(z)F 2 (z) so the 
exponents 7 r in F(z) are the sums of the exponents in F\(z) and in F 2 (z). This 
implies that K — K\ + K%. These different approaches to the asymptotics are 
reconciled by the classical but non-trivial fact that 

b(K x + K 2 ) = b{K x ) + b{K 2 ) 

for any two plane convex sets. This example may be used to construct counterex- 
amples to various conjectures. 

Example 9 Let us consider the entire function 

n 

F n (z) = n~ l J2 exp (e 2nir/n z) . 

r=l 

According to Corollary [7] 

N n (E) = b n E/2n + 0(l) 
as E — > oo, where b n —> 2tt as n — > oo. On the other hand 

lim FJz) = F(z) := — exp (e ld z) d9. 

n^oo ' 2lX h V ' 

Since F(z) is rotationally invariant and entire with -F(O) = 1 it must be identically 
equal to 1, and so it has no zeros at all. This establishes that even for this class of 
entire functions, the asymptotic distribution of the zeros does not vary continuously 
with the function. 

When combined with Corollary [?|, the following theorem proves that the eigenvalues 
of L move off to infinity as the eigenvalues of A coalesce. The singular behaviour 
of the spectrum in the limit 5 — » is in accordance with the behaviour for 5 = 
described in the Introduction. 

Theorem 10 Let A = aI + 5B, where a^O, the eigenvalues of B are all distinct, 
and 5 6 C is sufficiently small. Let L be defined in the usual manner with dim(C/) = 
p. Then there exist distinct constants a r such that 

lr = p (/3 - a ) f a + ar S + 0(5 2 ) (6) 

and there exists k > such that 

b(K) = k5 + 0{5 2 ). (7) 
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Proof The eigenvalues of A are a r = a + Sb r for 1 < r < n where b r are the distinct 
eigenvalues of B. The formula for the 7 r in Section 2 yields (|6]) immediately. This 
implies (|7]) with 

Q 

k = E |cr r — o r _ \/2ir. 

r=l 

In the remainder of this section we show how to apply Corollary [7] to certain elliptic 
systems with variable coefficients. We assume that 

Lf(x) = A(x)f'(x) 

in L 2 ((a,/3), C n ), where A(x) is invertible and diagonalizable for each x G (a,/?). 
In our next theorem we assume boundary conditions of the form f(a) = u and 
(f(/3),v) = 0, where m,d6 C n are both non-zero. 

The meaning of the word 'generically' below will be explained during the proof. 

Theorem 11 Let b(x) be the length of the boundary of K(x), where K(x) is the 
convex hull of the eigenvalues of A(x)~ l . If also A(-) is piecewise constant then 
generically one has 

N(E) = ^- f b(x)dx + 0(l) (8) 

Z7T J a 

as E — >• oo. 

Proof Following the notation of Section 2, equation (|3|) becomes 

F(z) = (V m e zD -V- 1 .V m . 1 e zD ^V~ 1 _ 1 ...V 1 e zD W 1 - 1 u,v) 

= (9) 

r=l 

Denoting the eigenvalues of A s by {a^t}^!, D s is a diagonal matrix with entries 
{a~j (a s — a s _i)}™ =1 . Each 7^ in (|9]) is of the form 

m 

^ = E a s~t(s)(«s - (10) 
s=l 

where t(-) is a function from {1, m} to {1, n}. The parameter r is a relabelling 
of the set of all such functions, so R < n m , with equality unless two sums of the 
form (0) happen to be equal. If K is the convex hull of the 7 r then our generic 
assumption is that 5 r 7^ for all vertices 7 r of K, and that no 7 r lie within any 
edge of K. 

By Corollary [7| we have to prove that 

b{K) = [ b(x)dx. (11) 

J a 

The right hand side equals 

m 

EK^Q(« S - a s-i) 

s=l 
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where A(x) = A s and K(x) = K s if a s -\ < x < a s as in Section 2. The identity 
(|TT|) follows by combining the facts that 

b{tK x ) = tb{K{) (12) 
b{K 1 +K 2 ) = b(Kx)+b{K 2 ) (13) 

for any convex sets K±, K 2 and any t > 0, with the identity 

m 

8=1 



Conjecture 12 We conjecture that Theorem |TT] remains generically valid for 
piecewise continuous coefficients A(x), x G [a, 0\. 

Note 13 The leading coefficient in @ depends only on the symbol of L. However, 
the truth of the theorem depends upon the genericity assumption for the following 
reason. Let L, L be two operators with the same symbol but different boundary 
conditions. Since the coefficients 5 r in the expansion ([§) may vanish for different 
values of r for the two operators, it may happen that K ^ K in an obvious notation. 
The leading coefficient in the spectral counting function @ will then usually be 
different. 

This phenomenon is illustrated by the following theorem. We make the same as- 
sumptions as in Theorem [TT], except that we now allow general boundary conditions 
of the form (|l|). The following theorem may be regarded as a version of Weyl's 
formula for the operator associated with the matrix-valued symbol 

where x G [oc,f3] and £ G R. Once again the meaning of the word 'generically' 
will be explained during the proof. We warn the reader that the asymptotic form 



obtained here is different from that of Theorem [Tl|; this is possible because the 
boundary conditions of that theorem are all non-generic in the sense used in the 
present theorem. 

Theorem 14 Under the above assumptions one generically has 

E 



N(E) = — ( b(x) dx + 0{ 

2ll Ja 



where 



b( X )=2Y: 



\a r (x)\ 1 



r=l 



and {a r (x)}" =1 are the eigenvalues of A(x) for each x G (a,/3). 
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Proof A more detailed understanding of (^) may be obtained by considering 
F( Zl ,...,z m ) = det(5 + TF TO e^ D ™y-^y m _ 1 e^- lD ™- 1 ^ 1 ...y 1 e^y 1 - 1 



R I m \ 

5^5 r exp X^^TV ■ ( 14 ) 

r=l Vs=l / 



The coefficient 7^ can be determined by putting 

f w if s = t 

2 = s 

otherwise 



to obtain 



R 

^ 5 r exp (u>777) = det (S + TV t e wDt V^ 1 

r=l 

= det (Vf 1 ) det (sV t + 7V t e wA 

By expanding the final determinant we deduce that each 777 is the sum of some 
subset of the numbers a^j(a t — a t -i) where 1 < j < n and {a tj }" =1 are the 
eigenvalues of A t . Since F(z) = F(z, z) we also have 

m 

8=1 

We say that we are in the generic case if all possible 5 r are non-zero, and also no 
7 r lies within an edge of the convex hull K of { , y r }^ =1 . 



We now apply Corollary 7 exactly as in the proof of Theorem |TT] but with different 
choices for K s . Generically each K s is the convex hull of all numbers of the form 

V s,J = a sJ 

je.r 

where J ranges over all subsets of {1, ...,n}. Hence 

K s = J S) i + ... + J S) „ 
where J SJ - is the interval with end-points 0,a~j. It follows by ([13|) that 



b(K s )=j2b(i SJ )=2Y: 

3=1 3=1 



a s,j 



-1 



This yields the statement of the theorem. 

We had expected that the above theorem would be valid for periodic boundary 
conditions, but the following example shows that this need not be the case. 
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Example 15 Let a = 0, (3 = 2, n 
for 1 < x < 2. Suppose also that 



2, A(x) = Ax for < x < 1 and A(x) = A 2 



ax,i 
a 2 i 



Ai = Vi 

for i — 1,2, where all are non-zero. Finally put 

n - u * 





Vr 1 



-i 



for z = 1, 2, where = a x , and Vi — u, 2i . 

Periodic boundary conditions correspond to the choice S = I and T = —I and 
lead to the formula 



F(z) = det(l-V 2 e D2Z V 2 - l .Vxe DlZ Vr 1 
= cdet(x -e D2Z Xe DlZ ) 

where c = det(V2) det(V A 1 _1 ) ^ and X = V 2 ~ 1 Vi is invertible. Hence 

[ ^(l-e^^) X li2 (l-e^^) " 



(15) 



If all X it j are non-zero we deduce that the possible values of 7 r are 0, U\ + u 2 , U\ + 
t> 2 , ui + u 2 , Vx + v 2 and Ux + u 2 + v i + t>2. Depending on the positions of these 
points K may have several shapes, and hence b(K) may have several values. One 
possibility is 

b(K) = 2|m 2 + ui| + 2|v 2 + wi| 

which is quite different from the formula obtained in Theorem H. 

In the very special case A 2 = —Ax every solution of Lf = zf is periodic and 
Spec(L) = C. But in this case u 2 + ux = v 2 + v\ = and X = I, so flOp yields 
F(z) = for all z. 



5 Basis Problems 

The fact that one may be able to determine the eigenvalues and eigenf unctions of 
an operator does not imply that the eigenfunctions form a basis in the relevant 
Banach space, 0, 0. If this fails then the positions of the eigenfunctions may 
be very unstable with respect to small perturbations of the operator, and the 
significance of the spectrum becomes moot. Boulton [||, || has investigated a closely 
related problem for second order elliptic systems in the language of pseudospectral 
theory, || [7| |[ [L5|, [16], pi}] , which amounts to estimating the resolvent norms of 



the operators concerned. 
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In this section we show that such problems do indeed occur in the context of this 
paper. For some operators of the type we consider the eigenfunctions form a basis, 
but for many others, and we believe most, they do not. We do not attempt a 
complete analysis, but just discuss the simplest example, of a first order system 
with two components. The operator H = L 2 provides a second order system with 
the same eigenvectors as L and hence the same basis problems. 

Example 16 Let n = 2 and O^iGR and let L be defined by 

(Lf)i(x) := tif[(x) 
(Lf) 2 (x) := -tif 2 {x) 

for all / G W 1,2 ((0, ir), C 2 ), subject to the boundary conditions /i(0) = /2(0) 
and /i(vr) = f 2 {^). A direct calculation shows that Spec(L) = tZ and that the 
eigenfunctions form a complete orthonormal set. Thus L is self-adjoint. 

The behaviour of the following slightly more general example is quite different. 
Example 17 Let n — 2, s > 0, t > and u := s + it. Let L be defined by 

(Lf)i(x) := uf[(x) 
(Lf) 2 (x) := uf 2 (x) 

for all / G jy 1,2 ((0, 7r), C 2 ), subject to the boundary conditions /i(0) = f 2 {fS) and 
/i(y) = ^(tt)- A direct calculation shows that A is an eigenvalue if 

e 7r\(i/u-i/u) = 1 

or equivalently if 

X n = (s 2 + t 2 )n/t 
for some nfZ. The corresponding eigenfunction is 

/n,i(x) = e xX -' u 



Lemma 18 Let S be a compact subset of C with zero Lebesgue measure and sup- 
pose that C\S has two components, one containing and the other unbounded. 
Then the linear span of the functions {z n } n ez is uniformly dense in C(S). 



Proof By Hartogs- Rosenthal lemma, |T(| p 47], C(S) = R(S) where the latter is 



the uniform closure in C(S) of the space of rational functions which do not have 
a pole on S. Each such rational function may be written as a linear combination 
of functions z m , (z — a)~ n where m, n > and o £ S. If a is in the unbounded 
component of C\S then (z — a)~ n may be uniformly approximated by polynomials 
in z by Runge's theorem, p 28]. If o is in the bounded component then 
(z — c)~ n may be uniformly approximated by polynomials in z~ x by using inversion 
and Runge's theorem. Putting these facts together completes the proof. 
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Theorem 19 The set of eigenf unctions {f n }nez is complete in the sense that its 
linear span is dense in L 2 ((0, tt), C 2 ). 

Proof Define w : [0, tt] — > C by 

w(x) = f 1A (x) = 

where t = s/t. Let S = T UT where T = w([0, tt]). Then S is a closed curve in C 
surrounding the origin and satisfies the conditions of Lemma [1^. 

It is sufficient to prove that if G C([0, tt], C 2 ) satisfies 0(0) = 4>{n) = then 
may be uniformly approximated by finite linear combinations of {/ n }nez- Given 
such a we define ip G C(S) by 

ij)(w(x)) = <pi(x) if x G [0, 7r] 
■0(w(a;)) = 02(s) if a; G [0, 7r]. 

Given e > there exists an approximation 

\\lp(z) - $rZ r \\oc < £ 

r=-N 



by Lemma [TBI Putting z — w(x) where x G [0, tx\ we obtain 



N 

1101 - XI ^r/r,l||oo < B. 
r=-N 



Putting z = w(x) where x G [0, ir] we obtain 

AT 

||02 - X <Wr,2||oo < B. 
r=-N 

Combining these we obtain the required estimate 

N 

110 - H ^r/r||oo < 2£. 
r=-JV 

The fact that the eigenvalues of L are real does not imply that it is similar to a 
self-adjoint operator. 

Theorem 20 Let P n be the spectral projection associated with the eigenvalue X n 
of L. Then \\P n \\ diverges at an exponential rate as n — > oo. The eigenf unctions 
f n therefore cannot constitute a basis of L 2 ((0, 7r), C 2 ). 

Proof We have 

p , _ ((f>,9n)fn 
\Jni 9n) 
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where g n is the appropriate eigenfunction of L*, and hence 

up || II 9n || 1 1 fn 1 1 

\{fn,9n) \ ' 

A direct calculation shows that the eigenfunction g n is given by 

g ntl (x) = ue' xX ^ 
g n ,2{x) = -ue~ xX "/ u . 

Evaluating the relevant integrals we find that 

{fn, 9n) = -2nit 

II f II 2 2-Ksn/t 

sn 

II M 2 / 2 i .2\ 

sn 

as n — > +oo, with a similar formula as n — > — oo. This implies the statement about 
||P n ||. The final statement is a consequence of the fact that if the eigenf unctions 
form a basis then ||P n || must be a bounded sequence, H, IS], [II . 



Note 21 If G L 2 and n = P n then it is possible that 

0~E<^ ( 16 ) 

neZ 

in the sense of some Abel-type summation scheme, |TJ], [TJ, |T^]. In the particular 
case in which <p{x) = (1, 0) for all x G (0, 7r) one finds that ||0 n || - * oo exponentially 
fast as n — > ±oo. It is not clear that the convergence of flip] ) using a summation 
scheme would have much numerical significance, because of the high instability of 
the spectrum under small perturbations of the operator. 



6 Numerical Range Problems 

The properties which we have discussed so far are similarity invariants. In other 
words if L has the property and T is a bounded invertible operator then TLT~ X 
also has the property. In this section we discuss properties which depend upon the 
particular norm chosen out of a similarity class. These are important because the 
norm is often given by physical considerations, and even if it is not, one might not 
know whether some better equivalent norm exists or how to find it. 

Let Hf = —Af" in L 2 ((0, C n ) where A is a bounded invertible n x n matrix. 
We impose any linear boundary condition at x — (3 and a boundary condition at 
x = of the form /(0) G U, /'(0) G V, where U, V are linear subspaces of C n . 

The following theorem, which extends results of Boulton, ||, [|, establishes that 
Kato's theory of sectorial forms, one of the main tools in non- self-adjoint semigroup 
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theory, cannot be applied to such operators even if the eigenvalues all lie in a half 
plane {z G C : Re (z) > k} for some k, as one might expect if the eigenvalues of 
A all have positive real parts. 

Theorem 22 If there exist c G U and d G V such that [Ad, c) 7^ then the 
numerical range of H equals the entire complex plane. 

Proof Let <p G C°°(R) satisfy <f>[x) = 1 if x < 1/3 and 0(x) = if x > 2/3. Let 
also i/)6R and n G Z + , and define 

/ n (z) = (c + e^d [(x + l/n) 2/3 - (l/n) 2/3 ] ) 0(x) 

so that / G Dom(i7). An easy calculation shows that 



lim ||/„| 



(c + e # dx 2/3 ) 0(x) 



dx 



which we denote by fci 7^ 0. An integration by parts establishes that 

(Hf n , f n ) = QUn) + B n 



where 



Q(fn) 



(Afn, fn) ^ 

(c + e^dx 2/3 ) <j>'{x) + § e^dx~ 1/3 0(x) 



dx 



as n — > 00, and 
We deduce that 



B n = (A/;(0),/ n (0)) = je^^c). 



(/m /n) ^1 

as n — > 00. By varying *0 G R and using the fact that the numerical range is a 
convex set, || Theorem 6.1], we deduce that it must equal C. 

Example 23 We consider the operator H defined on L 2 ((0, n), C 2 ) by 



(Hf)x 
(Hfh 



2ia fit 
e Jl 



-Ha t'l 



-e 

subject to the four boundary conditions /i(0) = 0, /^(O) = and 

cos(0)/ 1 (vr) + sin(0)/ 2 (7r) = 
-sm(0)#(7r)+cos(0)£(7r) = 0. 



To avoid redundancy we assume that — 7r/2 < 6 < tt/2. 
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In the particular case 6 = the two components of H are independent and its 
eigenvalues are n 2 e 2lOL where n = 1,2,... and also m 2 e~ 2ta where m = 0,1,2,... 
Moreover H is normal and its eigenfunctions form a complete orthonormal set in 
L 2 . 



If 6 7^ then we may apply Theorem |22| with c = (sin(#), — cos(#)) and d = 
(cos(6 l ), sin(6 1 )). We obtain 

(Ad, c) = i sin(20) sin(2a) 

which is generically non-zero. We deduce that the numerical range of H equals 
C. This fact cannot, however, be discovered simply by finding the eigenvalues. 
If 6 7^ then is not an eigenvalue of H. The eigenfunction associated to the 
eigenvalue A = z 2 must be of the form 



fi (x) = Ci sinh (. 
f 2 (x) = c 2 cosh( 



xze 



xze ia 



if it is to satisfy the first two boundary conditions. The existence of a non-zero 
eigenfunction satisfying the other two boundary conditions then forces F(z) = 
where 

F(z) := e ia cos(8) 2 smh(irze- ia )smh(iTze ia ) 

+ e~ ia sm(6) 2 cosh(7rze" 4a ) cosh(7rze iC *). 

This is of the canonical form 



F(z) = E^ e 



r=l 



where 71 = 27rcos(a), 72 = 27rzsin(a), 73 = — 27rcos(a) and 74 = — 27risin(a). 
The points 7 r are the vertices of a rhombus, and the directions of the four outward 
pointing normals are ±ie ±ia . For \a\ < 7r/4 this implies that all except finitely 
many of the eigenvalues A = z 2 lie in the half plane C = {z : Re (z) < 0}. We 
conjecture that all the eigenvalues of H satisfy Re (A) < 0. 



7 Semigroup Properties 

In spite of the above, there are cases in which one can prove that the operator H 
is the generator of a contraction semigroup on L 2 ((a,j3), C n ) provided this space 
is given a new equivalent norm. The following two theorems describe the abstract 
situation. 

Theorem 24 Let Z be a closed operator on the Hilbert space TC and let A be a 
bounded invertible operator. If ^ Spec(Z) and AZ is accretive then —AZ is 
the generator of a contraction semigroup and hence Spec(AZ) C C + , the set of 
complex numbers with non-negative real parts. 
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Proof The operator AZ is closed and invertible, so there exists e > such that 
z ^ Spec(AZ) for all \z\ < e. The theorem now follows by applying || Theorem 
2.25]. 

Theorem 25 Let Z be a closed operator on the Hilbert space TC and let A be a 
bounded invertible operator. If Z + XI is invertible for some X > and AZ, A are 
both accretive, then —AZ is the generator of a contraction semigroup and hence 
Spec(AZ) C C+. 

Proof By the hypotheses the operator A(Z-\-XI) is invertible and accretive. Hence 



—A[Z + XI) generates a contraction semigroup by Theorem gj. Since A is a 
bounded perturbation, — AZ generates a strongly continuous semigroup. This 
must be a contraction semigroup by || Theorem 2.27]. 

We now return to a more concrete context. Let Z = — d 2 / da; 2 act in L 2 ((a, j3), C n ) 
subject to the boundary conditions f(a) G U, f{(3) G U, /'(a) G V , f'(/3) G V 
where dim([7) + dim(V) = n. 

Theorem 26 IfUnV ^ {0} then Spec(Z) = C. IfUnV = {0} i/ien Z zs szmz/ar 
to a self-adjoint operator with non-negative discrete spectrum. 

Proof Let ^ c £ ([/ U ^)" L - Then for any z G C the function 

(/(*) = ce rzx 

is orthogonal to the range of Z — 2 2 /, so z 2 G Spec(Z). 

If, on the other hand, U H V = {0} then we may write Z as the (non-orthogonal) 
direct sum of its restrictions to L 2 ((a, (3), U) and L 2 ((a, P),V). In the first we are 
imposing Neumann boundary conditions, and in the second Dirichlet boundary 
conditions, so the spectrum is as stated. Another way to express the same idea is 
to give C n an inner product which makes U and V orthogonal, so that Z becomes 
self-adjoint with respect to this inner product. The change from one inner product 
to an equivalent one amounts to the same as a similarity transformation. 

Note 27 In the above theorem G Spec(Z) unless U = {0} and V = C n , the full 
Dirichlet case. This is why we need both of Theorems |24| and 



Note 28 If the angle between U and V is very small then the condition number of 
the similarity transformation is large and it will be difficult to distinguish between 
the two alternative conclusions of the above theorem. In other words Z will have 
bad pseudospectral properties. 

The following provides a partial converse to Theorem We have not been able 
to find a similar result if the boundary conditions are different at the two ends of 
the interval. 
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Theorem 29 Let Hf = —Af" in L 2 ((a, (3), C n ) where A is a bounded invertible 
n x n matrix. We impose boundary conditions f(a) G U , f{(3) G U , /'(a) G V , 
f{0) G V" where U f) V = {0} and [/ + V = C n . T/ien if A is accretive and 
AV _L U, the operator —H generates a contraction semigroup in L 2 ((a, j3), C"). 

Proof We can apply Theorem with Zf = —f" subject to the above boundary 
conditions provided we show that H is accretive. If / G Dom(H) then integration 
by parts yields 

(HfJ}= f{Af,f)dx. 

J a 

This has non-negative real part since A is accretive. 

The condition AV _L U depends upon the particular inner product used in C n and 
is not satisfied generically. On the other hand the condition AV D U — {0} does 
hold generically. The second condition implies the first for a suitable choice of the 
inner product on C ra , but one then needs to check whether A is accretive for this 
new inner product. Assuming this does happen, if the condition number relating 
the two inner products is k then one only gets 



with respect to the standard inner product for all t > 0. This may be of value if k 
is fairly small but it is of little computational use if k is sufficiently large. 

Example 30 One can prove that L generates a one-parameter group in the follow 
situation. Let Lf = Af subject to the quasi-periodic boundary conditions f{(3) = 
Sf(a), where A is a diagonal matrix with real, non-zero eigenvalues and S is 
invertible. It appears that L generates a periodic flow on [a,/3], but with different 
flow rates in the different components and mixing at the end points. One can also 
regard this as a flow in a generally non-trivial vector bundle over the circle. We 
restrict attention to the simple case A = I, but note that Boulton's ideas allow 
one to treat the case in which the eigenvalues of A are real and rationally related, 

Under the above assumptions we have 

F(z) = det {S - e z{l3 ~ a) l} . 

Putting a = e z ^~ a ^ and finding the eigenvalues of S we obtain n generally distinct 
solutions a = a r . We conclude that the eigenvalues of L lie on one of n distinct 
lines parallel to the y-axis. 

Our next theorem implies that it is only possible to prove that e Lt is a strongly 
continuous one-parameter group using numerical range ideas if L is already skew- 
adjoint. In spite of the failure of this method, the subsequent theorem shows that 
L does generate a one-parameter group whenever the eigenvalues of A are all real. 



19 



Theorem 31 The numerical range of L is only contained in a strip of the form 

{z e C : ci < Re 0) < c 2 } 

if L is already skew adjoint, in other words if S is unitary. 

Proof For general functions in Dom(L) we have 

Re(LfJ) = ||/(/3)|| 2 -||/(«)|| 2 

= ((S*S-I)f(a)J(a)). 

If 5 is not unitary then this may be arbitrarily large (positive or negative or both) 
for functions in Dom(L) of unit norm; the proof is similar to that of Theorem |22|. 
This prevents the numerical range from lying in a strip of the stated kind. 

Theorem 32 For all operators L defined as in Example L is the generator of 
a strongly continuous one-parameter group acting on L 2 (a, f3),C n ) . 

Proof We define the bounded invertible operator W : L? — > L? by 

Wf(x) = S(x)f(x) 

where S(x) is any smooth function with values in the invertible n x n matrices 
satisfying 5(a) = / and S(0) = 5. Now let V be the space of functions in W 1 ' 2 
such that f{a) = f(0) and define M : V -> L 2 by 

Mf = W^LWf. 

A direct calculation shows that 

Mf(x) = f'(x)+Y{x)f{x) 

where Y is a bounded smooth matrix- valued function on [a, f3\. Writing this is the 
more abstract form M = Mq + 7 we see that M is a bounded perturbation of Mo, 
which generates a one-parameter group of isometries on L 2 . Hence M, and then 
L, generate one-parameter groups on L 2 . 

8 Inverse Spectral Theory 

In 0, |J, Boulton showed that a second order elliptic system H may have real 
spectrum without being similar to a self-adjoint operator. In this section we inves- 
tigate the extent to which one can reconstruct a first order elliptic system from its 
spectrum. We start with a partial positive result and then give an example which 
shows that in general one cannot go any further. 

We continue in our earlier framework including the generic assumption of Section 3. 
The first condition of the following theorem is probably unnecessary. We conjecture 
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that in general the multiplicity of any eigenvalue of L equals the order of the 
corresponding zero of F(z). Note that we have already proved that all zeros of 
F(z) of large enough modulus are simple, so the condition only concerns a finite 
number of smaller eigenvalues. 

Theorem 33 Assume that the zeros of F(z) are all simple. Then the spectrum of 
L determines the function F(z) up to an exponential factor, and hence the set of 
constants 7 r , 1 < r < R, up to a common additive factor. 

Proof Although we do not use this in the proof, we note that by Theorems [5] and 
^| the asymptotics of the spectrum of L determine the length and direction of each 
edge of K. Since K is a convex polygon it can only be reconstructed in one way 
up to translations in C. 

The statement of the theorem is a consequence of the fact that F{z) is an entire 
function of order 1, and can therefore be written in the form 



n=l 



F(z) = z m e hz J] { (l - y) eZ/Zn \ ( 17 ) 



for some h G C, where z n are the zeros of F(z). See [12] , p. 199]. The fact that 



F(z) determines all of the 7 n 1 < r < R, up to the additive factor h is elementary. 

Note 34 One may break the product in QT7D into sporadic terms together with 
terms associated with z n which are close to one of the straight lines already de- 
scribed. The products associated with each line are asymptotically similar to 
certain formulae involving gamma functions, namely 



e (7+^ r ( a + 1) °° 



r(a - z + l) „=i l \ a + n 



z 



e z/(a+n) 



where 7 is Euler's constant and 



£i (a + n)n 

This follows immediately from Whittaker and Watson, E3, p 236]. 



In certain cases one can go further than Theorem |33|. If dim(VF) = 1 then we 
showed in Section 2 that 7 r = a) /oJT for 1 < r < R where a r are the eigenvalues 
of A. If in addition we know the centre of K and the value of /3 — a, then we can 
evaluate all of the eigenvalues of A. The following example shows that in general 
one cannot go any further. 

Example 35 Given k > let u = s + it lie on the circle C = {u : s 2 + t 2 = kt) and 
consider again Example [17]. We exclude the case s = t = for obvious reasons; 
the case s — 0, t — y/k is covered by Example |l^. The two vertices of K are 
7± = 7r(u ± ik) where v = s/(s 2 + 1 2 ) = s/kt. We thus see that K is the same, up 
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to a translation, for all of these operators as u varies on C. By Examples [16] and 



17| the spectrum of L is kZ for all choices of u on C. 

In order to see that L cannot be reconstructed from its spectrum up to a similarity 
transformation, we need to prove the following theorem. 

Theorem 36 No two of the operators of Example W\ are similar as u varies in C. 



Proof We use the notation of Theorem 120. If two such operators L and L de- 



termined by u, u G C were similar then it would follow that ||P n ||/||P„|| would be 
bounded above and below by positive constants uniformly with respect to n. 



Now the proof of Theorem |20] implies that 



~7vknv L. I 

I p I 1 * l " 
± n 



2n\t\knv 



as n — > oo. We deduce that if L and L are similar then v = v. This implies that 
u = u. 



Acknowledgments I would like to thank J M Anderson, L S Boulton, K M Ball, 
D A Edwards, G Rozenblum, S G Scott and R F Streater for valuable discussions. 



References 

[1] Agranovich M S: Elliptic Boundary Problems. Encycl. Math. Sci., p 1-140, 
Springer (1997). 

[2] Aslanyan A, Davies E B: Spectral Instability for some Schrodinger operators. 
Numer. Math. 85 (2000) 525-552. 

[3] Bottcher A: Pseudospectra and singular values of large convolution operators. 
J. Int. Eqns. Appl. 6 (1994) 267-301. 

[4] Boulton L S: Ph D thesis, King's College, London, 2001. 

[5] Boulton L S: Spectral Behaviour of a simple non-self-adjoint operator. 
Preprint, 2001. 

[6] Davies E B: One-Parameter Semigroups. Academic Press, London, 1980. 



22 



[7] Davies E B: Pseudospectra of differential operators. J. Oper. Theory 43 (2000) 
243-262. 

[8] Davies E B: Semi-classical states for non-self-adjoint Schrodinger operators. 
Commun. Math. Phys. 200 (1999) 35-41. 

[9] Davies E B: Wild spectral behaviour of anharmonic operators. Bull. London 
Math. Soc. 32 (2000) 432-438. 

[10] Gamelin T W: Uniform Algebras. Prentice-Hall Inc., Englewood Cliffs, NJ, 
1969. 

[11] Gohberg I C, Krein M G: Transl. Amer. Math. Soc. 18 (1969) 309-316. 

[12] Hille E: Analytic Function Theory, vol. 2. Blaisdell Publ. Co., Waltham, Mass., 
1962. 

[13] Lesch M. Tolksdorf J: On the determinant of one-dimensional elliptic boundary 
value problems. Commun. Math. Phys. 193 (1998) 643-660. 

[14] Lidskii V B: The Fourier series expansion in terms of the principal functions 
of a non-selfadjoint operator, (in Russian) Mat. Sb. 57 (1962) 137-150. 

[15] Reddy S C: Pseudospectra of Wiener-Hopf integral operators and constant 
coefficient differential operators. J. Int. Eqns. and Applic. 5 (1993) 369-403. 

[16] Reichel L, Trefethen L N: Eigenvalues and pseudoeigenvalues of Toeplitz ma- 
trices. Linear Alg. and its Applic. 162-164 (1992) 153-185. 

[17] Scott S G, Torres F: Relative zeta functions and relative determinant bundles 
on manifolds with boundary: (I) dimension one. Preprint, 2000. 

[18] Scott S G: Zeta determinants on manifolds with boundary. Preprint, 2000. 

[19] Shkalikov A A: Elliptic equations in Hilbert space and associated spectral 
methods. J. Soviet. Math. 51 (1990) 2399-2467. 

[20] Streater R F: Stability of a hot Smolochowski fluid. Open Systems and Infor- 
mation Dynamics. 7 (2000) 1-9. 

[21] Trefethen L N: 1997 Pseudospectra of linear operators. SIAM Review. 39 
(1997) 383-406. 

[22] Whittaker E T, Watson G N: A Course of Modern Analysis. Cambridge Univ. 
Press, Cambridge, 1927. 



Department of Mathematics 
King's College 



23 



Strand 

London WC2R 2LS 
England 



e-mail: E.Brian.Davies@kcl.ac.uk 



24 



